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ABSTRACT 


Hydrographic data has traditionally been used only to determine dynamical 
features, specifically dynamic topography and geostrophic currents relative to an assumed 
level of no motion. However, the assumption of geostrophic balance has been questioned 
because of the neglect of ageostrophic effects. Measuring vorticity and vertical motion 
has always been technically demanding and costly. Due to the importance of these 
parameters new techniques have been developed, among than a diagnostic technique 
known as the C-vector Method, which uses the wind stress and the geostrophic balance 
as forcing functions to infer these quantities. 

Using this technique, the three dimensional vorticity field was computed for the 
California current system during the Coastal Transition Zone (CTZ) program (17-26 
March, 1987). Hydrographic (CTD) and wind data sets were used, assuming a "quasi- 
geostrophic" system. Also computed was the vertical vorticity of the C-vector, known 
as the T function, which has a direct relationship with the vertical velocity, making it 
possible to infer upward or downward motion in coastal waters. 

A satisfying correlation was found between surface temperature fields, satellite 
imagery, and vertical motions inferred from theT function. The C-vector has therefore 
been shown to be a very reliable method of diagnosing vorticity and vertical circulation. 
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I. INTRODUCTION 


A. PURPOSE 

Direct measurement of ocean currents with sufficient spacial and temporal coverage 
to be useful is technically demanding and costly. In shallow water, current meters may 
be ‘ "pended from an anchored ship in order to measure currents at several depths 
simultaneously. However, this procedure only gives information about the currents at a 
single geographic location. Also, an anchored ship does not remain stationary, but rather 
moves in relation to the anchor. This is a source of error for the measured current that 
increases with the length of the anchor line. Another, more practical method is to 
suspend current meters from a moored buoy. A number of such buoys moored in a 
pattern in the ocean will provide information about the three dimensional distribution of 
currents as a function of time. However, because of the expense, difficulties at sea, and 
the complicated nature of the currents, it is not possible to obtain observations over as 
much of the ocean as we would like. Thus, the geostrophic method is the primary 
method to determine currents. 

The geostrophic method calculates the currents indirectly from the density 
distribution in the ocean; it is easier to measure temperature and salinity than to measure 
currents directly. In fact, most of our knowledge of the ocean circulation below the 
surface has been obtained in this way. The geostrophic method is also useful in strong 
currents for which it is difficult to moor recording current meters. This method still has 
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some disadvantages, such as die lade of vertical motion information (Pond and Pickard, 
1983). 

This work deals with the development of a new sy stem for diagnosing the three- 
dimensional ageostrophic circulation in the Coastal Transition Zone (CTZ) using 
hydrographic and wind data sets. In this system, die flow contains two parts: geostrophic 
circulation and ageostrophic circulation. First, the geostrophic currents are obtained 
directly from the hydrographic data using traditional methods. Then, the ageostrophic 
circulation is assumed to be forced by die geostrophic currents (a phenomenon called 
geostrophic forcing) and by the surface wind stress. 

R THE CALIFORNIA CURRENT SYSTEM 
1. Water Masses 

There are four water masses in the California current system. Each has its 
own characteristic temperature, salinity, dissolved oxygen, and nutrients. 

a. Pacific Subarctic Water 

Following Simpson (1984), this is formed from the Oyashio and Kuroshio 
extension and carried eastward to the North America continent as part of the Subarctic 
Current and West Wind Drift. This flow collides with the continent near 48N. Part of 
it separates and is diverted south into the California current with water that is relatively 
low in temperature and salinity, but high in oxygen and phosphate. 

Although this water mass mixes with die California Current and increases 
its temperature as it moves southward, it gives the offshore California Current its low 
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salinity and high dissolved oxygen near die surface. 

b. North Pacific Central Water 

Formed in the central Pacific between depths of 100 to 800m, the North 
Pacific Central is the least saline of the central water masses of die Pacific Ocean. This 
warm, salty water is relatively low in dissolved oxygen and nutrients, and alters die 
California System from the west. 

c. Upwelled Waters 

Found within 50km of the coast, this water is relatively cold, salty, and 
high in nutrient content and oxygen. Upwelled water is the result of the Honan transport 
of surface waters away from the coast at times of equatorward winds along the coast. 

d Equatorial Pacific Water 

Formed in the eastern tropical Pacific below a strong thermocline with 
waters of high temperature, salinity and nutrients, but low dissolved oxygen. Equatorial 
Pacific water enters die California Undercurrent from the south (Simpson, 1984). 

2. Current Systran 

The California Current System consists of the California Current and the 
California Undercurrent. 

The California Currents lie between 0 and 200m depth. Offshore, the flow 
takes place equatorward during die whole year. Near the coast there is a fall-winter 
reversal, creating the California Countercurrent (Simpson, 1984). 
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TABLE 1.1. RELATIVE COMPARISON OF WATER MASSES OF THE 
CALIFORNIA CURRENT SYSTEM 
(L = Low; H = High) 

(Simpson, 1984) 



T 

S 

Qz 

Nutrients 


Surface Water Masses 
(0-200m) 





a 

Pacific Subarctic 

L 

L 

H 

H 

b. 

N. Pacific Central 

H 

H 

L 

L 

c. 

Coastal Upwelled 

L 

H 

L 

H 


Subsurface Water Masses 
(200-500m) 





d. 

Equatorial Pacific 

H 

H 

L 

H 


The California Undercurrent, with its core at about 300m depth, originates in 
the eastern equatorial Pacific and runs poleward throughout the year. This poleward flow 
has beat observed at different levels in different seasons (Huyer et al., 1989). 

In the winter season there is a broad northward surface current, named the 
Davidson Current. It is seen between October and March, with the strongest mean flow 
in January or February at a speed of20-30 cm/sec. North of Pt. Conception the Davidson 
Current is wider than the shelf but the strongest flow is over the inner shelf against the 
coast. The current is forced by the southeasterly winds in winter, which cause 
down welling at the coast due to the Honan transport. This combines with increased 
coastal runoff during the rainy season and sets up an offshore density gradient that causes 
the northward geostrophic current (Huyer et al., 1989). 
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Also, there are inshore surface currents during the upwelling season in 
southern and central California These currents seem to be inshore of die upwelling jet 
that flows south. The poleward flow is evident most of die year in die Southern 
California Bight. In that area it appears to form the inshore side of a large cyclonic eddy 
that fills the Bight. Further north, the poleward surface current can appear whenever the 
equatorward winds relax. It is much narrower here, 10-2Okm compared to up to 100km 
in the Bight. Velocities reach about 30 cnVs. A similar current is off central Oregon in 
mid and late summer It is believed that these countercurrents are driven by an along¬ 
shore pressure gradient that balances the equatorward wind stress. That is why they are 
most intense when the winds relax (Huyer et al., 1989). 

3. Wind Forcing 

a. Hind Driven Circulation 

The wind driven circulation affects only the upper few hundred meters. 
Therefore, it is considered predominately a horizontal circulation in contrast with the 
thermohaline circulation. 

b. Coastal Upwelling 

Coastal upwelling occurs along the entire west coast of the North 
American continent, and is primarily due to Ekman transport of surface waters offshore 
forced by the alongshore equatorward winds . 

This upwelling is forced by atmospheric circulation around the North 
Pacific High which varies seasonally in both strength and position. The High migrates 
from its most southerly position at 28°N, 130 °W in February, to its most northerly 
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TABLE 1.2. SUMMARY OF THE POLEWARD FLOW NEAR NORTHERN 
AND SOUTHERN BOUNDARIES OF THE U.S. WEST COAST 

(Hickey, 1989) 


WASHINGTON SHELF/SLOPE VS. OREGON SHELF/SLOPE 


Washington Shelf Washington Slope Oregon ShelfrSlope 


Spring/Summer 


Equatonvard at all Poleward (below 
depths (spring). 100m). 
Considerable vertical Sub-surface 
shear and flow maximum at shel 
reversals over upper break. 

20m Vertical shear 


to 400m) due to 
e strong equatorward 

at shelf winds here vice 
Washington. 




maximum 


Poleward 

Sub-surface 


iiliiili. it.. 


winter due to 
Columbia River 
plume). 


Poleward but less 
* so than Washington, 

at --400m Decreases from mid 
to outer shelf. 
Equatorward 
vertical shear. 

Some flow reversals 
with depth 


SOUTHERN CALIFORNIA BIGHT SHELF VS. SLOPE 



Spring/Summer 


Fall/Winter 


Shelf 

Eastern Santa 
Monica Basin Slope 

Weston Santa 
Monica Basin Slope 

Equatorward except 

Poleward from 40m 

Poleward 

at 30m and 60m 
during August (for 
narrow shelf in 
southern part of 
bight). 

ot 100m everywhere 
except off San 
Onofre in spring at 
40m 

Sub-surface 

maximum 

Sub-surface 

maximum 

Vertical shear 
maximum like 
Washington slope. 


Poleward 

Equatorward (after 


strongest during 
these seasons. 


everywhere except 
off San Qnofre at 
40m in early fall. 
Sub-surface 


October at 40m and 
100m) 


maximum 




































Atmospheric Pressure 
(mb) 

JANUARY 



Atmospheric Pressure 
I (mb) 

I JULY 


10125 


1005 




m 


H 

'102$ mb 


j^?a^ mn Term mean a,mos P he ™ 


sea level for January and 
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Figure 1.2. Mean wind stress 1854-1972. Contours are constant magnitudes of 0.5, 
1.0, and 1.5 dynes cm 2 ; in shaded areas the wind stress exceeds 1.0 (Nelson, 1977). 
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position at 38°N, 150°W in August, and bade again; the maximum pressure increases 
from wintt. (1020mb) to summer (1025mb) as shown in Figure 1.1. The atmospheric 
pressure over the extreme western United States is much higher in winter (-1020mb), 
when the surface is cold, than it is in summer (~1015mb), when die surface is warm. 
This thermal effect on pressure is particularly marked in California's central valley, where 
the Thermal Low reaches 1005mb in July. The offshore pressure gradient along the west 
coast is strongest in summer. The strength of the pressure gradient across coastal waters 
varies with position along the coast, with season, and obviously with day to day variations 
in the weather. 

As shown in Figure 1.2, the region of maximum equatorward wind stress 
migrates northward from 25°N in January to 39°N in July, as expected from the pressure 
distribution shown in Figure 1.1. The strongest equatorward winds are observed in July 
off northern California, where the center of the High is closest to the Thermal Low. 
South of about 40°N, winds are equatorward all year round. North of this latitude, the 
wind stress is northward (e.g., unfavorable for coastal upwelling) in winter. 

The alongshore wind stress, T, causes a net offshore transport in the 
surface Ekman layer, 11f, which in turn equals the amount of water that upwells along the 
coast. The Coriolis parameter, £ increases from 0.6 x lfrV at 24°N to 1.1 x lO'V at 
about 5(fN; thus favorable winds of the same magnitude will result in nearly twice as 
much upwelling off southern Baja California as off northern Washington. Short-term 
fluctuations with periods of several days in the wind stress can cause upwelling, even 
when die mean wind stress for a given month is not favorable. 


9 





Open ocean upwelling can occur away from die coastal boundary if there 
is a positive curl in the wind stress field, i.e., if there is net divergence in the Ekman 
transport field For example, if offshore winds are more strongly southward than coastal 
winds, there is stronger seaward Ekman transport offshore than near the coast, and there 
must be upwelling away from the coast as well as along the coast (Figure 1.3). Gradients 
in the actual wind field are likely to be very much stronger near shore, particularly near 
large headlands such as Cape Mendocino, and along steep coastal mountain ranges 
(Huyer, 1983). 
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D. DATA 


Hydrographic (LTD) data used in this study was collected from a cruise off central 
California during 17-26 March 1987 as part of the ONR Coastal Transition Zone (CTZ) 
Program. The CTZ program studied die spring summer transition in the California coast. 
The spring transition is defined as the event that causes die passage from the typical 
winter regime with level isopycnals and no mean vertical shear to the typical 
spring/summer regime with sloping isopycnals and a mean southward surface current and 
a strong mean vertical shear (Huyer et al., 1979). It was observed to occur in mid-March 
during the CTZ in 1987. 

A. Coastal Transition Zone Program 

The CTZ program was designed to understand the dynamics and implications of 
structures known as "filaments." Filaments are cold features in surface temperature, 
typically less than 100 km wide but hundreds of kilometers long which extend offshore 
from near die coast. Typically, they are accompanied by high concentrations of 
chlorophyll. 

Easton boundary currents have been characterized as eddy-rich and strongly time- 
dependant with spatially structured variability in the flow field. Bernstein et al., (1977) 
demonstrated that the complex structures visible in satellite infrared (IR) imagery could 
be matched to comparable variability measured throughout traditional in-situ means. 
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Satellite images revealed long tongues or filaments of cold, high-chlorophyll water 
extending from the coastal zone out to 200 km from shore. These were associated with: 

1. strong, narrow seaward currents (surface drifters suggested upwelled waters 
directed offshore), 

2. continuation of a strong shoreline coastal jet flowing generally southward along 

the upwelling front; and 

3. pairs of oppositely rotating eddies. 

Four cruises were made during February/June 1987. Buoys, ships, satellites, 
aircrafts and several instruments were used to collect diverse variables and data (Brink 
and Cowles, 1991). 

& Hydrographic Data 

The hydrographic data was collected from a cruise off central California during 17- 
26 March 1987 as part of a study of the CTZ project, and was reported by Jessen et al. 
(1989). The study area was from Pt Reyes, California to Pt. St. George, Oregon from the 
coast to 150 km offshore. 

1. Planned Sampling Grid 

The planned sample grid shown in Figure 2.1 consists of two shoreline 
sections 150 and 90 km offshore, a discontinuous shoreline section 60 km offshore, and 
four across-shore sections. The across-shore sections divided the sampling grid into three 
separate modules (A,B, and Q. 
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2. Actual Sampling Grid 


The actual sampling grid differed from the original due to weather and time 
constraints. A total of 96 CTD casts (Figure 2.2) to 500 m depth and 55 XBT drops 
(Figure 2.3) were made during this cruise. 

Additionally, maps of hourly averaged wind vectors, surface temperature, 
salinity, and dynamic height relative to 500 db are presented in Figue 2.4-2 7, 
respectively. 

3. Data Acquisition and Calibration 

Hydrographic data was acquired using a Neil Brown Mark III-B CTD and 
Sippican T-4 XBTs. A General Oceanic Rosette Sampler was attached to the CTD and 
equipped with twelve 5 liter Niskin bottles for in-situ water sampling. The CTD sampling 
rate was 32 Hz, but the acquisition software employed a latch filter which limited the 
number of data points collected during each cast to 4,308. On the 500 m casts this 
resulted in the acquisition of 8 or 9 data points per meter of water CTD data was 
acquired only on the downcast with a winch speed of approximately 30 rrVmin to 150 m 
then 60 m/min to 500m The data was acquired using an HP200 computer and stored on 
3.5 inch diskettes. Upon return, the data was then transferred to a nine track tape and 
processed on an IBM 3033 mainframe computer. In addition to the CTD and XBT data, 
an underway data acquisition loop recorded 30-second averages of sea surface temperature 
and salinity, sea surface skin temperature, wind speed and direction, air temperature, and 
dew point temperature. The sensors used to acquire this data included Seabird 
temperature and conductivity sensors for the sea surface temperature and salinity, a 
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Figure 2.1. Planned CTD 

(Jessen, et al., 1989). 
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Figure 2.2. Actual CTD station numbers and locations for Cruise CIZ1 during 
March 17-26,1987 (lessen, et al., 1989). 
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Figure 2.4. Hourly average of wind speed and direction measured at 10m height 
on March 87. (Jessen, et al., 1989) 
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Rosemount 100-ohm platinum resistance thermometer for the sea surface skin temperature, 
a Young anemometer for the wind speed and direction, and a General Eastern dewpoint 
sensor for the air and dewpoint temperatures. The underway data was acquired chi an 
HP9816 computer and recorded on 3.5 inch diskettes. Like the CTD data, the underway 
data was transferred to a nine track tape upon return and processed chi die IBM 
mainframe The temperature, conductivity, and pressure sensors on the CTD and the 
temperature and conductivity sensors of the underway sampling system were calibrated 
shortly before the cruise. The pressure calibration was carried out using a Chandler 
Engineering dead weight tester as a standard. The temperature calibration was done using 
a Seabird temperature sensor as a standard. The conductivity calibration was carried out 
using a Guildline Model 8400 Autosal as a standard. 

A total of 61 water samples were taken at 8 CTD stations for post-anise 
calibration. The pressure, conductivity and temperature were noted as each sample was 
taken. These numbers, after applying the pre-cruise calibration coefficients, were used 
to calculate salinity and the results compared with the water sample salinities calculated 
using the Guildline Model 8400 Autosal in the laboratory. 

4. Data Processing 

After the raw CTD data was transferred to the IBM 3033 mainframe computer 
at the Naval Postgraduate School, temperature and conductivity corrections were applied 
to produce profiles of corrected pressure, temperature, and conductivity. Salinity was 
calculated from these corrected values according to the algorithm of Lewis and Perkin 
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(1981). Severe spiking due to system malfunctions was eliminated from die salinity 
signal with a search for vertical salinity gradients great® than 1.0 PSU/m Points that 
were determined to be bad were replaced by the values computed from the linear 
interpolation. Time lag spikes were eliminated by discarding salinity data in regions 
where the vertical temperature gradient exceeded 0.2° C/m and replacing the discarded 
data with linearly interpolated values. Finally, the data was averaged within 1 m intervals 
and visually examined for any remaining outlier missed during processing. If found, these 
points were replaced with linearly interpolated values (Jessen et al., 1989). 
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ID. C-VECTOR METHOD 


As mentioned before, ocean currents are die result of die combined effects of 
thermohaline motion and of wind-driven motion. In both cases, the induced motion 
usually continues far beyond the place where it was initiated. It should be appreciated 
that while the wind primarily causes horizontal motion, there may also be vertical motion 
such as the upwelling and down welling both in the open ocean and near the coast. 
Similarly, while the thermohaline effects primarily give rising and sinking motions 
(vertical), then by the continuity, horizontal motion will be induced (inflow or outflow). 

Hydrographic data of coastal oceans were traditionally used to determine dynamical 
features in only two fields: dynamic topography relative to a certain depth, and 
geostrophic currents by assuming a certain level of no motion. Some information is lost 
when using these methods. 

In reality, the ocean flow is not in absolute geostrophic balance. It contains two 
types of circulations, the geostrophic current and the ageostrophic current. The former 
is obtained from the hydrographic data sets with the usual technique (geostrophic 
balance), and the ageostrophic current is induced by the presence of the geostrophic flow 
and surface winds. Hence, the three-dimen si onal circulation which is both geostrophic 
and ageostrophic can also be calculated by the hydrographic and surface wind data sets 
(Chu, 1992). 
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In regions where the temperature and salinity gradients are strong, die usual 
procedure to calculate only geostrophic currents gives possible serious errors. 

To determine the three dimensional flow, a new "tool" is presented since neither die 
geostrophic current nor dynamical height can supply enough information about the three- 
dimensional flow field for the coastal currents. 

A. BASIC THEORIES 


In Meteorology many papers were published for diagnosing the vertical motions in 
the atmosphere. One of them is the Q-vector method. Based on the "Quasigeostrophic 
Theory," the quasigeostrophic CO equation was first presented by Hiassen (1962) and 
Hoskins (1978). To fill gaps in Hoskins's theory, Xu (1992) proposed a C-vector concept 
for the atmospheric mesoscale moist firontogenesis. The C-vector concept has been 
adapted to physical oceanography (Chu, 1992) for diagnosing three-dimensional 
ageostrophic flow fields. The geostrophic currents (Ug,Vg) are computed as. 


= 3 Li£ = -ito 

dz p 0 dy dy 


(3.1a) 


& v g _ _ g dp _ db 
dz p 0 dx dx 


(3.1b) 


where ft, is the characteristic sea water density, p is the density deviation from pb, R, + 
p is the observed density field, and b = -gp^pb is the buoyancy force. 





However, these geostrophic currents are not the complete currents. The basic 


equations for the complete three dimensional currents are: 

(3.2a) 

{ Tt*V’ V)V * f{U - U 9 )S *W <3.2b) 

{j£+V-V)b+N 2 w=0 (3.2c) 


du dv dw 
dx Uy dz 


=0 


(3.2d) 


where V; v = (u,v,w), f is the Coriolis parameter, N is the Brunt-Vaisala frequency, and 
(Y*,y>) ; turbulent momentum flux caused by the surface wind stress (x K ,'f). The three- 
dimensional flow (V) can be broken down into geostrophic (\Q and ageostrophic (V. 
) components. 


Y *\ 




(3.3) 


V = V g +V a (3.4) 


Thai, after the decomposition, the equations (3.2a)-(3.2d) (basic equations) should 
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have die following structure: 


3(v a ) = 9t(v g , P, y x , yy) (3.5) 

where 3 , 9? are differential operators. 

Thus, we can use this equation (3.5) with the hydrographic data (Y p) and surface 
wind data, to calculate the three-dimensional ageostrophic circulation, if the real forms 
of these operators are acquired. 


R QUASI-GEOSTROPHIC SYSTEM 
In the quasi-geostrophic system, we have 


a 

at 


+ v • v 


at 


+ v g -V 


(3.6) 


The basic equations for the coastal water (3.2a)- (3.2d) on an f-plane become 


fV a = 


dY* 

dz 


- (-Y- +V 

' at 


V) u a 


(3.7a) 


fu a 


dyy _ ,_a 
dz at 



• V) v g 


(3.7b) 


N 2 w ’ 


= - [JL+v 


at 


v g *)b 


(3.7c) 


du a dv a atv a 
dx dy dz 


= 0 


(3.7d) 
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Three-Dimensional 
Ageostrophic Circulation 


Figure 3.1. Schematic diagram for the computation of ageostrophic circulation(Qiu, 

1992). 


Applying the cross-derivative technique cm the ageostrophic pseudo-vorticity 
equations, 

~jz (f2v * ] + '£p {N * w * ) =2Ci (3 * 8a) 

-^(f 2 u a ) ~^(N z w a ) - 2 C 2 (3.8b) 

-A(f 2 v a ) --A(f 2 u a ) = 2C 3 (3.8c) 

Rearranging these equations, we get the components of the C-Vector, 


r W*dV g dV'dU'^f e>Y* 
1 ' dy dz dy dz ; 2 dz 2 


(3.9a) 
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(3.9b) 


C 2 




dU g dV 
dz dx 




+ 1H2X1 

2 dz 2 


du g dv g dv g du g f d , dY* A dyy, 

2^ l “5F 


(3.9c) 


Qearly we can see that die C-Vector, C = (Q.C^Q), is computed from die 
geostrophic currents (l^, Vg) (obtained from the hydrographic data set), and the wind data 


(v,n 


Figure 3.1 shows the diagnostic model for die three-dimensional ageostrophic 
circulation. 


C PHYSICAL MEANING OF C-VECTOR 
1. Ageostrophic Vortex Line 

The following nondimensional form for flows with f (coriolis parameter) and 
N (Brunt-Vaisala frequency) as constants, gives the mathematical connection between the 
C-vector and the ageostrophic pseudo vorticity equations (3.7a-c). 

VxV a =2R 0 C (3.10) 

V-? a =0 (3.11) 


Ro s U/(fL) is the Rossby number and H/L = W/U = fTN is chosen for die scaling where 
H is the diaractistic vertical length, L is the characteristic horizontal length, W is the 
characteristic vertical velocity, and U is the characteristic horizontal velocity. 

Since the ageostrophic vorticity is proportional to the C-vector, a C-vector 






streamline can be viewed as an ageostrophic vortex line (Figure 3.2). 



Figure 3.2. Ageostrophic circulation induced by C-Vector in three-dimensional 
space (Xu, 1992). 

2. Vertical Motion Induced by the Vertical Vortidty of C-Vector 
The vertical velocity equation is obtained from 

V 2 !/ a = -2J? 0 £-VxC (3.12) 

which indicates that the vertical velocity is caused by the vertical vorticity of the C- 
Vector (Figure 3.3.). 

3. Nondivergence of C-Vector 

To have nondivergence, the horizontal divergence (convergence) should be 
balanced by the vertical convergence (divergence), the summation of differentiations of 
(3.9a) with respect to x, of (3.9b) with respect to y, and of (3.9c) with respect to z, leads 
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Figure 33. The vertical vorticity W induced fay horizontal rotation (vorticity) of C- 
Vector (Xu, 1992). 

to 


V-c=o 


(3.13) 


which demonstrates that the C-vector is nondivergent. Additionally, (3.8c) shows that the 
vertical vorticity of the ageostrophic flow can be deduced from the vertical component 
of the C-vector (Figure 3.4) 




V a 

to 

Figure 3.4. The vertical vortidty of ageostrophic flow is inferred from the 
vertical component of die C-Vector (Xu, 1992). 

D. 4* Function 

As defined before in equation 3.12 there is a direct relation between the vertical 
velocity and the vertical vorticity of the C-vector, which it can be rewritten as. 


V 2 w a = -2R 0 ¥ 


(3.14) 


where 


¥=£*VxC 


(3.15) 


This 4* function is very important in the determination of the vertical motion of the 
flow because both should have the same sign i.e., if 4 / is positive it will have upward 
motion, a downward motion when 4* was negative. 
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E. C -Vector 


The C -Vector defined in die following relation as: 

C* s C-VECTOR + 3-DIMENSIONAL GEOSTROPHIC VORnCITY 

CT =(0,(^0) 

It is also derived from the C-vector to give: 


, _ Ci _ dv g 
x ~ f2 dz 


(3.16a) 


c s.£k + 2 flifa 
y f 2 2 dy 


(3.16b) 


z f 2 2 dx 


dv„ du g 
dy 


(3.16c) 


This means that this vector C = (Q, G,, CJ, is another and very useful way to 
determine the 3-dimensional vorticity for the total flow. 
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IV. THREE DIMENSIONAL PSEUDO-VORTIOTY FIELD 

The averaged data was extracted using die IBM mainframe computer at the Naval 
Postgraduate School using die fortran program "EXTRACT" written by Paul Jessen. This 
program extracts data from a file cm die mass storage system and writes it to a file on the 
user's reader disk for easy access. 

A. Interpolation Using The Weighting Function 

Since the extracted data was not regularly distributed and evenly spaced, it was 
necessary to put it in a form that makes future calculations easy. Hence, an interpolation 
was performed, creating a new distribution of die data on a new grid (figure 4.1). This 
interpolation was done using the so-called "weighting function” technique. 

This technique was developed in order to achieve an optimal interpolation on data 
that has samples (CTD casts as an example) that are relatively coarsely distributed This 
analysis technique was developed by Doswell (1977) and Maddox (1980), who also based 
it in a Barnes (1964/1973) scheme (filtering properties of the objective analysis). As 
explained by Tintore et al. (1991), die Barnes scheme has two analyses. The first one, 
f b defines the macroscale, a low pass filter which smooths the observed data, the second 
one, f 2 , another low pass filter with a cut-off near the smaller resolvable scale, then the 
mesoscale analysis is defined as the normalized difference between these two low-pass 
filters (a band-pass filter). Finally, the total field is obtained with the sum of these two 
scales (macroscale and mesoscale). At the same time, each analysis consists of two steps 






explained using the following formulas. 


a. First Step 






_ / 3»1 


n *1 


-j) 


and 


w n (i,j) =exp [ 


-d n (i,j) 2 
4 C 


(4.1) 


(4.2) 


Where f°(i j) are grid values initially computed from the total N of data observed; f are 
data observed, W^(ij) is a weight function; d„(ij) is the separation distance between the 
grid point (ij) and then* station; and C is an analysis parameter. 

b. Second Step "The First Guess" 


j) +-2^- (4.3) 

E 
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(4.4) 


A f B =f a -fZ 

iy n -(i,j). e xJ : d ^ )i ] 0< G <1 (4.5) 

where £° is the first guess evaluated at each station; and G is another analysis parameter. 

The analysis parameters C aid G are chosen either as a result of previous 
knowledge of structures in die field cm- based upon theoretical constraints such as 
characteristic scale and deformation radius. 

R Grid Presentation 

The extracted data was subdivided into two areas, in order to facilitate the 
interpolation work. 

Area 1 coincides with section A in the planned grid Figure 2.1. It is formed 
from 27 stations (CTD casts). After performing the interpolation. Area 1 has 49 grid 
points (x=7,y=7), as shown in Figure 4.1. In this area, the X-axis is zonal (positive 
eastward), and Y-axis is meridial (positive northward). 

Area 2 is formed with 56 stations corresponding to sections B and C of 
planned grid Figure 2.1. After the interpolation. Area 2 has 225 grid points (x=9,y=25). 
For this area the Y-axis coincides with the shoreline direction, and X-axis with the cross¬ 
shore. The Z axis is in the vertical direction. The maximum depth in this study is 500m. 


36 




The vertical resolution is 10m, giving a total of 50 grid points in the Z direction In 
summary. Area 1 is formed by (x=7,y=7,z=50), and area 2 by (x=9,y=25,z=50). 


C Geos trophic Velocity 

The method to compute the geostrophic velocity was performed on the interpolated 
data (grid point). This method was written based on the thermal wind equation where 
the horizontal density differences between the stations are used to compute the baroclinic 
shear (assuming a level of no motion of 500 m). Then the shear is vertically integrated 
to give the geostrophic velocity relative to 500 m 

To explain this method better we will assume we want to calculate the geostrophic 
velocity between stations A and B, with density profiles p a (z) and p,(z) respectively, 
following these steps: 

a Compute an average value of the Coriolis parameter, £ for the two stations. 

f.+f* 


f= 




(4.6) 


b. For k = 1 to NPTS-1, compute a horizontal gradient at level k + 1/2 


dp _ [Pfe(fc)+P fc (k+1)~P«(k)-p a (k+1)] 
dx 2Ax 


(4.7) 


where NPTS is the number of grid points in z direction and Ax is the distance between 
the two stations. 

c. Compute the baroclinic shear, from the thermal wind equation 
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4r=-_2_4£ 

oz fp 0 ox 


(4.8) 


at level k + 1/2, where g = 9.81 m/i and ft, = 1026 kg/n 1. 

d. Assuming a level of no motion at the lowest depth (500m) set V(NPTS) = 0, and 
then integrate die baroclinic shear upwards in z to compute V at level k. 

V(k) =V(k+l) + lz(lc)-z(lc+l)) (4.9) 

Now die question should be, how do we know if both the interpolated method and 
the geostrophic computation (based in the thermal wind equation) are correct? The 
answer is shown by the correspondence found between plots made using die program 
PLOTSEC, also written by Paul lessen, for die geostrophic velocity computed by the 
GEOVEL program (with hydrographic data), and die plots generated by die subroutine 
(CONTOURS for die geostrophic velocity computed using die thermal wind method 
(interpolated data). Figures 4.2.a-b, 4.3.a-b, 4.4.a-b, 4.5.a-b, 4.6.a-b, show these plots for 
Area 1, and Figures 4 7.a-b and 4.8.a-b for Area 2. 

It is important to note that the boundaries of die plots do not coincide exacdy since 
one is for real station positions, and the other is the grid position. Another important 
note is that die grid points are evenly spaced and perfectly aligned, but the stations are 
neither evenly spaced nor perfectly aligned. Despite these differences, they show a good 
correlation. 


38 




D. Wind Stress 


The bulk equation that relates the stress! to the wind speed u has the form 

t = Cpp u 2 (4.10) 

where Q, is die drag coefficient (dimensionless) and depends on the roughness of the 
surface. 

The drag coefficient C D for die ocean increases with the wind speed, and it was 
calculated for the California Current System from Large and Pond (1981) using 

C fl =l.l4xl0 -3 (4.11) 

for u < 10 m/s and 

10’ 3 <^=0.49+0.06 5u (4.12) 

for u > 10 m/s. 

Noting that die wind stress only affects the upper layers of the ocean, the so-called mixed- 
layer depth (ML) for this study was based on the data and 65m was assumed. 
Furthermore, it was assumed that the wind stress is maximum at the ocean surface, 
decaying linearly to zero at the base of die ML. 

E. C-Vector 

We can now compute the C-vector by applying die equations presented in Chapter 
III and generate a series of plotted fields for presentation and analysis in the n ext chapter. 
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Figure 4.1. Final grid distribution of data after interpolation, using the weighting 
function technique. 
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Figure 4.2.b. Geostrophic velocity from Area 1, Cross-Section y=l, computed using 
THERMAL WIND method with interpolated data. 
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Figure 4.3.b. Geostrophic velocity from Cross-Section v=4, Area 1, computed using 
THERMAL WIND method with interpolated data. 
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Figure 4.4.b. Geostrophic velocity from Area 1, Cross-Section y=7, computed using 
THERMAL WIND method with interpolated data. 
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Figure 4.5.a. Geostrophic velocity from Area 1, Stations 88-54, computed using 
GEOVEL with hydrographic data. 
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Figure 4.5.b Geostropbic velocity from Area 1, Cross-Section x=l, computed using 
THERMAL WIND method with interpolated data. 
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Figure 4.6.a. Geostrophic velocity from Area 1, Stations 92-50, computed using 
GEOVEL with hydrographic data. 
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Figure 4.6.b. Geostrophic velocity from Area 1, Cross-Section x=6, computed using 
THERMAL WIND method with interpolated data. 
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Figure 4.7.a. Geostrophic velocity from Area 2, Stations 56-27, computed using 
GEOVEL with hydrographic data. 
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Figure 4.7.b. Geostrophic velocity from Area 2, Cross-Section x=l, computed using 
the THERMAL WIND method with interpolated data. 
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Figure 4.8.a Geostrophic velocity from Area 2, Stations 50-01, computed using 
GEOVEL with hydrographic data 
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Figure 4.8.b. Geostrophic velocity from Area 2, Cross-Section x=6, computed using 
the THERMAL WIND method with interpolated data. 
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V. RESULTS 


A. PRESENTATION 

1. Pseudo-vorticity Field for the Total Flow 

Q corresponds to the total pseudo-vorticity in the x-direction. There are 
seven Q cross-sections in Area-1 (Figures 5.1-5.3) and nine Q cross-sections in the Area- 
2 (Figures 5.7-5.9). Q corresponds to the total pseudo-vorticity in the y-direction. There 
are seven Q cross-sections in Area-1 (Figures 5.4-5.6) and twenty five cross-sections 
in the Area-2 (Figures 5.10-5.15). The depth of all of these cross-sections is 500m, with 
10m vertical resolution. The principal characteristics for the analysis are: solid lines for 
positive and zero pseudo-vorticity (i.e. clock-wise rotation), dashed lines for negative 
pseudo-vorticity (i.e. counter clock-wise rotation). 

2. VP Function 

The plots of the Tfunction are presented in Figures 5.16-5.18 for Area-1 and 
in Figures 5.19-5.21 for Area-2. In these horizontal cross-sections the solid lines 
represent upward motion and the dashed lines downward motion. Zero lines are also 
represented by solid lines. 
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E ANALYSIS 

1. Area-1 

In all die plots presented for this area significant values of magnitude for Q 
Cy and ¥ Function were found from the surface to a depth of 150m, but die level of 
greatest magnitude seems to be at a depth of 100m . 

In Figure 5.17 the ¥ function is plotted at three different depths, 60, 80 and 
100m. At the 100m depth, we can see upward motion (solid lines) and downward motion 
(dashed lines) at sectors listed in Table 5.1. 


| TABLE 5.1 VERTICAL MOTION FOR ¥ (Area 1) 

X-Coordinates 

Y-Coordinates 

Vertical Motion 

x=l-4 

y=l-3 

Upward 

x=6-7 

y=4-7 

Upward 

x=l-4 

y=6-7 

Upward 

x=l-4 

y=4-5 

Downward 

x=4-7 

y=l-3 

Downward 

x=5 

y=6 

Downward 


This analysis is very similar in that it shows upward motion. Indeed, even 
in Figure 5.18 the same characteristics are present at the 120m level. Deeper than 120m, 
die features change, and their intensity decreases and almost disappears at the 160m level. 
To complete the analysis, each component of the pseudo-vordcity (QQ) will be 
calculated with the ¥ function. 
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a. C x and T Function 

In Figure 5.1 three vertical cross-sections are presented for Q at x=l-3. The 
pseudo-vorticity field is clearly depicted, positive clock-wise (solid lines) and negative 
counter clock-wise (dashed lines). The cores of rotation (both positive and negative) are 
found mainly around the 100m depth. The same depictions are found in Figures 5.2 and 
5.3 for C x and in Figures 5 4-5.6 for Q,. To analyze and compare the pseudo-vorticity 
fields with the 'F function, we choose the 100m depth as a base level. 

In Figure 5.17 the western boundary that corresponds to x=l as constant and 
moves along the y axis from south to north (y=l-7). Compared with Figure 5.1 at x=l 
(also at 100m depth), this corresponds to the vertical cross-section of Q As we can see 
at y=l, a core of positive vorticity is found. Going to the left at y=4 a core of negative 
vorticity is found. Thus, we expect upward motion in between, going even farther to the 
left near y=7, a weaker center of positive vorticity is found, giving downward motion to 
the left of the center of negative and to the right of positive vorticity. Finally, for this 
section we find upward motion to the left of the core of positive vorticity, which agrees 
with the T function at this level 

In Figure 5.1 at x=2, only two centers are present, negative to the right at y=3 
and positive to the left at y=6. As depicted by the arrows upward motion at y=l-3 and 
y=6-7, and downward motion at the center at y=3-6 are found. At x=3 a similar analysis 
for x=l is done, where upward motion is at y=l-3 and y=6-7, downward motion at y=3-6. 

In Figure 5.3 the analysis is very similar, but interesting features are found 
at x=5; at y=l-3 and y=6-7 downward motion is found and upward motion at the center 






y=3-6. In sections x=6 and x=7 both of the single positive vorticity cores are at die 
center y=4. Then, upward motion is found to die north from y=4-7 and negative to the 
south from y=4-7. This analysis agrees with the'F' function, as anticipated. 

b. C y and Function 

In Figure 5.4 three vertical cross-sections are presented for at y=l-3. The 
pseudo-vorticity field are clearly depicted, positive clock-wise (solid lines) and negative 
counter clock-wise (dashed lines), the cores of rotation (both positive and negative) are 
again found mainly around the 100m depth. The same is true for Figures 5.5 for y=3-5 
and 5.6 for y=5-7. To analyze and compare the pseudo-vorticity fields with theT 
function we choose the 100m depth as a base level. 

In Figure 5.17 the south boundary that corresponds to y=l and moves along 
die x axis from west to east (x=l-7). Compared with Figure 5.4 at y=l, this corresponds 
to the vertical cross-section of (^, as we can see at x=5, a core of positive vorticity is 
found, and at x=l a core of negative vorticity. We expect upward motion in between 
these centers, and downward motion at x=5-7. Now, at y=2, almost the same situation 
as in y=l is present giving to upward motion at x=l-4 and downward at x=4-7. At y=3 
only two centers are present, positive at x=5 and negative at x=2. Then, as depicted by 
the arrows we found upward motion in between both colters at x=3-5 and downward to 
at x=l-3 and x=5-7. 

In Figure 5.5 and 5.6 the analysis is very similar. As expected, this analysis 
agrees with the HP function and could be repeated at any level and with any cross section. 
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2. Area-2 


The technique to analyze this area is similar to that applied in Area* 1. All of 
the plots presented for this area for Q, C> and ^Function show activity from the surface 
until near 150m depth, but die level of greatest activity seems to be around the 80m 
depth. 

In Figure 5.19 the 'Ffunction is plotted at three different depths, 40, 60 and 
80m, where we can see upward motion as solid lines and downward motion depicted by 
dashed lines. Deeper than 140m (Figure 5.21), the magnitude is not great. Also, we can 
see that the greatest magnitude is found horizontally to the north and especially to the 
northwest, with very little to the south and almost nothing at center. 

3. Remarks 

It is interesting to note some important results on the levels of greatest activity 
in both areas where the cores of positive and negative vorticity and the areas of maximum 
upward or downward motion were found to affect the baroclinic zone, areas of strongest 
gradients on temperature, salinity, and density. This is especially true when compared 
with areas of upward and downward motion with the density fields. The tilt of the 
isopycnals suggests vertical motion at this level as shown in Figure 5.22-24, vertical 
cross-sections equivalents to y=l, y=4 and y=7, respectively, for Area-1 as an example. 

Another interesting point is that in Area-2 there was more activity at the north 
than at the south, and almost nothing at the center. The map of dynamic height (Figure 
2.7) shows the strongest features are at north and south and not too much in the center. 
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C CORRELATION 


1. Sea Surface Temperature 

Figure S.2S shows a map of surface temperature which overlapped die grid 
used for this study. In this map we can see a strong correlation between areas of 
upwelling and the V Pfunction at die lm depth. As an example, for Area-1 some filaments 
of cold water are present in different sectors such as x=l-3 and y=l-2, x=l-3 and y=6-7, 
x=5-7 and y=5-6, all of which are in agreement with the'F function at die 1 m depth as 
shown in Figure 5.26. For Area-2 the cold filaments are found at x=6-9 and y=18-22, 
x=2-9 and y=14-16, x=8-9 and y=5-9, also in agreement with the 4function at lm depth 
as shown in Figure 5.27. 

2. Satellite Imagery 

Data obtained from the Advanced Very High Resolution Radiometer (AVHRR 
cm - AVHRR/2) instruments carried on the TIROS-N/NOAA series of polar orbiting 
satellites yield the brightness temperature images for the study region. All AVHRR 
instruments measure emitted radiation in four wavelength bands, visible (0.6-0.7 pm), near 
infrared (0.7-1.1 pm), and thermal infrared (3.5-3.9 pm and 10.5-11.5 pm). NOAA-9 
AVHRR data was left as black body brightness temperatures from channel 4 (10.5-11.5 
pm) and were uncorrected for atmospheric contamination. The satellite data were first 
calibrated and navigated to earth coordinates. Each image was then co-registered to an 
identical coordinate map of the study region. 

Figure 5.28 also shows important information that correlates our analysis. 
This satellite photography corresponds to the 27 of March at 2200Z. Gear areas represent 
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upwelling and darker areas represent downwelling Inferred from die satellite images, an 
analysis of upward and downward motion in Area-1 is found at sectors is listed in Table 
5.2. 


TABLE SJ CORRELATION BETWEEN SATELLITE IMAGERY AND 4* 

(Area 1) 

X-Coordinates 

Y-Coordinates 

Circulation 

Correlation with T 

x=l-4 

y=l-2 

Upward 

Yes 1 

x=l-4 

y=i-2 

Upward 

Yes 

x=l-3 

y=5-7 

Upward 

Yes 

x=l-4 

y=2-5 

Upward 

Yes 


As we can see, all of these sectors have good correspondence with the^F 
function in Figure 5.26, although part of this area is obstructed by a patch of thick clouds. 

In the same way, Area-2 upward and downward motion was found in the 
sectors listed in Table 5.3. 
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TABLE 5J CORRELATION OF SATELLITE IMAGERY AND ¥ (Area 2) 


X-Coordinates 

Y-Coordinates 

Circulation 

Correlation with 4* 

x=T-7 

y=20-25 

Upward 

Yes 

x=7-9 

>'=20-22 

Upward 

Yes 

x=3-5 

y=17-19 

Upward 

Yes 

x=6-8 

y=16-17 

Upward 

Yes 

x=2-9 

y=l3-16 

Upward 

Yes 

x=6-9 

y=l 1-13 

Upward 

Yes 

x=7-8 

y=6-ll 

Upward 

Yes 

x=2-7 

y=4-6 

Upward 

Yes 

x=7-9 

y=4-5 

Upward 

No 

x=2-5 

y=2-4 

Upward 

Yes 

x=l-2 

y=l-2 

Upward 

Yes 

x=3-5 

y=i-2 

Upward 

Yes 

x=7-8 

y=l-2 

Upward 

Yes 

x=8-9 

y=23-25 

Downward 

Yes 

x=l-5 

y=16-20 

Downward 

Yes 

x=6-9 

y=17-20 

Downward 

Yes 

x=l-2 

y=8-12 

Downward 

Yes 

x=4-6 

y=5-9 

Downward 

Yes 

x=7-9 

y=5-6 

Downward 

Yes 

x=2-5 

y=4-5 

Downward 

Yes 

x=4-6 

y=3-4 

Downward 

Yes 

x=6-9 

y=3 

Downward 

Yes 

x=2-4 

y=l 

Downward 

Yes 

x=6 

y=l 

Downward 

Yes 


There exists a good correlation between areas of great variability and areas 
where there is not much variability. It is important to notice that during this time 
some features may move somewhat. In general, however, the analysis demonstrates 
consistency between the satellite information and theT function. 
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a CONCLUSIONS 

a) A three dimensional pseudo-vorticity field can be obtained from this new concept 
C-vector method using only CTD measurements (hydrographic data) and wind data. We 
can either directly visualize three dimensional circulations from the pseudo-vorticity field, 
or solve the equations (3.8a,b,c) to get three dimensional flow field. 

b) Also using this method we can infer vertical motions especially from the close 
relation obtained between the vertical velocity with the'F function. 

c) At this time, the "quasi-geostrophic" system assumed in this study is considered 
to be a valid approach to the problem, but it does not accurately represent conditions 
found in the real world. Development of an approach that considers conditions left 
unaddressed by the quasi-geostrophic assumption will enable even more realistic results 
to be achieved . 
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CTZ1 »- 1 

DENSITY ANOM (y) 10 KM 


54 53 52 51 50 



Figure 5.22. Vertical section of density anomaly (Jessen, 1989). Compare with y=l 











CTZ1 

DENSITY ANOM ( 7 ) 

83 82 81 80 79 



Figure 5.23. Vertical section of density anomaly (Jessen, 1989). Compare with v=4 
from Area-1. 










Figure 5.24. Vertical section of density anomaly (Jessen, 1989). Compare with y=7 
from Area-1. 
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Figure 5.25. Map of surface temperature Cruise CTZ1 (Jessen, 1989), overlapped 
by the grid used for this study. 
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Figure 5.27. 'F function at lm depth from Area-2. 
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Figure 5.28. Satellite (NOAA-9) AVHRR infrared image of brightness temperature 
(March-27-2200z). lighter shades indicate colder brightness temperatures. 
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